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Big bang nucleosynthesis in a modified gravity model of f{R) oc is investigated. The only 
free parameter of the model is a power-law index n. We find cosmological solutions in a parameter 
region of 1 < n < (4 + ■\/6)/5. We calculate abundances of '*He, D, ^He, ’^Li, and ®Li during big 
bang nucleosynthesis. We compare the results with the latest observational data. It is then found 
that the power-law index is constrained to be (n — 1) = (—0.86 ± 1.19) x 10“^ (95 % C.L.) mainly 
from observations of deuterium abundance as well as ^He abundance. 
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I. INTRODUCTION 

The present standard cosmological model is based 
on Einstein’s general relativity with the Friedmann- 
Lemaitre-Robertson-Walker metric for a homogeneous 
and isotropic universe. All elementary particles of the 
standard particle model and dark matter and dark en¬ 
ergy are taken into account in the cosmological model. 
The standard cosmological model has been supported 
by various kinds of astronomical observations. Obser¬ 
vations of light element abundances in old astronomical 
objects are, however, one of the most important premises 
of the standard cosmological model. Roughly speaking, 
the theoretical predictions of light element abundances 
are consistent with observational data. In modified grav¬ 
itational theories, cosmic expansion histories are different 
from that in the standard model, while in modified par¬ 
ticle theories additional effects of exotic particles operate 
in the early universe. As a result, primordial elemen¬ 
tal abundances in these models are different from those 
in the standard big bang nucleosynthesis (BBN) model. 
Therefore, we can limit any models which predict changes 
in abundances. 

The baryogenesis in a modified gravity model of 
f{R) oc i?", where R is the Ricci scalar and n is the 
power-law index, has been studied to explain the small 
baryon-to-photon number ratio of the Universe [l|. The 
authors derived a cosmological solution in which the scale 
factor of the universe scales as a{t) oc where t is the 
cosmic time and a is a real parameter. They argued that 
(4 - y6)/5 < n < 1 should be satisfied in order to re¬ 
alize a positive temperature of the universe. The BBN 
in the same model has also been analytically studied Q ■ 
They constrained the index to be 1 — n ;< 2 x 10“"^ by a 
comparison of an analytical estimation of ^He abundance 
and observational data. 

In this paper, we calculate BBN in the model of 
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f{R) cx i?" with a detailed nuclear reaction network 
code and show abundances of all light elements produced 
during BBN. In Ref. Q, only ^He abundance has been 
studied semianalytically. In this paper, however, it is 
found that observational constraints on the primordial 
D abundance can limit the modified gravity model more 
stringently than those on the ^He abundance. On the 
other hand, limits derived from observations of ^He, ^Li, 
and ®Li abundances are less stringent than those of D 
and "‘He. In addition, we point out that models of f{R) 
should describe the accelerated expansion of the present 
Universe. We find that the model used in the previous 
study P, @1 is excluded by this requirement, and we sug¬ 
gest a simple correction to the model. In this paper, we 
consider three models: (1) a new model which describes 
the accelerated emansion of the present Universe, (2) the 
previous model P, l2| which cannot describe the expan¬ 
sion, and (3) a corrected version of (2) which describes 
the expansion. Although the limit on the f{R) oc i?" 
model is corrected, our revised result supports the pre¬ 
vious conclusion that the consideration of BBN excludes 
parameter values of n largely different from unity 0. 

In Sec. n the modified gravity model is introduced, 
and equations for the cosmic evolution are derived. In 
Sec. mi our code for the BBN calculation is briefly ex¬ 
plained. In Sec. lYl observational constraints on the 
primordial light element abundances are described. In 
Sec. El a result of BBN is shown and interpreted. In 
Sec. ED this work is briefly summarized. 

II. COSMOLOGY OF f(R) oc R" GRAVITY 

In this section, formulas of the cosmology in the modi¬ 
fied gravity model are shown. First, we derive equations 
of motion. The action is given by 

5 "=^/ d‘^Xy/^f{R)-\-Smigf,u,(l}m), ( 1 ) 

where = SttG is defined, with G Newton’s constant, 

the metric tensor, g the determinant of the metric 
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tensor, and the action of the matter field (j)^ which 
takes into account radiation and matter in the Universe. 
The field equation for gravity is then derived by varying 
this action with respect to the metric tensor, 


The second model is given [T^ by 

/2(i?) = -(^)”. (11) 


fR^.u - If9^. - V^v,/' + g^,Df = (2) 


where f = df/dR is defined, and is the energy- 
momentum tensor for matter defined as 


T _ 2 6 gd^m) 


(3) 


Here is the Lagrangian density of matter, and it is 
related to S'm by 


Sm= J d^xV^C^. (4) 

We assume the spatially flat Friedmann-Lemaitre- 
Robertson-Walker metric, as supposed in the standard 
cosmological model, 

ds^ = dt^ — a{t)^ {dx^ + dy^ + dz^) . ( 5 ) 

For matter, on the other hand, we assume a perfect fluid 
described by a time-dependent energy density p{t) and 
pressure p{t), 

= diag (p, -p, -p, -p). (6) 

The 0-0 component of Eq. ([2) then becomes 

-5-f'-lf + 3-f"R = K^p. ( 7 ) 

a 2 a 

The i-i components, on the other hand, give 

+2^] f + 1/ - 2 “/"R- r- f'R = n^p. (8) 
\a J 2 a 

From Eqs. 0 and (Ill) , the following energy conservation 
holds, 


p-k3-(p-kp) = 0. (9) 

a 

Second, we constrain the model space to be studied in 
this paper. We assume two functional shapes for f{R). 
The first model is given 0,0 by 

MR) = (f) , (10) 

where H is a constant given by H with 

Mp = 1.22 X 10^® GeV the Planck mass. The power-law 
index n is the only free parameter, and n = 1 reduces 
to Einstein’s general relativity. This model has been an¬ 
alyzed in Refs. 00, and proper solutions of the scale 
factor exist only for (4—•\/6)/10 < a < 1/2 in this model. 


A formulation of this model is given using the same as¬ 
sumptions as those adopted in Refs. 0 , 0 , as follows. It 
is assumed that the matter part is predominantly con¬ 
tributed by the radiation with p = p/3. In this case, Eq. 
(ini) leads to a relation of p oc a~^. Here, we additionally 
constrain the model space by assuming the power-law 
solution of the scale factor, i.e., 

a(t)oct“. (12) 

Inserting Eqs. (ITT]) and (TT^ into Eqs. ([7]) and ([5]), we 
have two equations. 


. (2n + a — 3) 1 


2 (2a- 1) 




6a (2a — 1) 


At^ 


(13) 


and 


n (3a — 1) ^ n (n — 1) (2a — 2u-k 1) 1 

6(2a- 1) 3a (2a - 1) 2 


= K®p 


6a (2a- 1) 


We note that the Ricci scalar is given by 

6a (2a- 1) 


i?= -6 



/ • \ 2l 

d 

f a\ 

- + 

- 

a 

^aj 




(14) 


(15) 


It is found that a = u/2 must be satisfied in order to 
hold Eqs. m and (fTTI) for any time t. Then we assume 
a = n/2 in what follows. In this case, the energy density 
is related to the cosmic time as 

1 6a (2a — 1) 

^ [ A^ 

The pressure satisfies the relation p = p/3. Then we 
utilize the relation between the energy density and the 
cosmic temperature, 


— lOa^ -k 8a — 1 
2(2a- 1) 


(16) 


P=^9.T\ (17) 

where T is the temperature and g*{T) is the relativistic 
degrees of freedom for energy density. From Eqs. m 
and (113, the time-temperature relation of the universe 
is derived as 


T = 


15 y/" 1/4 


(18) 


where 


52a — 


(6a)"“ 


— lOa^ -k 8a — 1 
2 (2a- 1)^”^“ 


(19) 
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is defined. Since the cosmic temperature must be posi¬ 
tive, the parameter 52 a must be positive. Then a con¬ 
straint on a is derived, 


1 4-h Ve 

- < a < -—. 

2 - - 10 


( 20 ) 


because of nonreal fi{R) values for i? < 0 and n 1. 
The latter model is therefore excluded. 

The model fi{R) can, however, be corrected so that it 
is consistent with observational evidence of accelerated 
expansion of the present Universe. One simple correction 
is given by 


The Hubble expansion rate is given as a function of en¬ 
ergy density by inserting Eq. CH) into the equation 
H = aja = a/t. We thus obtain 


H = ag^ 


-l/(4a) 

06 


Mp 




( 21 ) 


Finally, the Hubble expansion rate can also be directly 
given as a function of temperature using Eq. dlB 


H = 


Hl/2 /4^3^^y/(4«) 


92. 


l/(4a) 


M, 


l/{2a) ( 

n \ 


15 


/C6 


( 22 ) 


/(i?) = sguiR) 


(24) 


This function reduces to fi{R) [Eq. (fTUl) ] for R>0 and 
f 2 {R) [Eq. pip ] for R < 0. This function also describes 
the general relativity (n = 1) in which f{R) = R/A is 
satisfied. It includes both solutions for a > 1/2 and 
a < 1/2, and allows a transition from i? > 0 to i? < 0 
in the cosmic evolution. Therefore, the model f{R) is 
consistent with the acceleration of the Universe. We then 
use this function in the following calculation. When this 
model is used, the time-temperature relation is given by 
Eq. (fT5P with the function g 2 a replaced by 


The special case of n = 1 (a = 1/2) corresponds to Ein¬ 
stein’s general relativity. We note that the allowed region 
[Eq. (j^ j is outside the region, (4 — •\/6)/10 < a < 1/2, 
in the /i (R) model 0,0. 

The two models [Eqs. (TTOP and dTlTl ] successfully de¬ 
scribe solutions of i? > 0 and R < 0, respectively. In 
the present setup, the left-hand sides of Eqs. m and (|5]) 
are proportional to the function /(i?). Since fi{R) oc i?" 
should be real for real number n, the scalar R must be 
non-negative. This requirement on R is important when 
we derive the constraint [Eq. (|20P ]. For the case of f 2 {R), 
on the other hand, the scalar R must be non-positive. 
Thus, in this model the power-law functions fi{^ and 
f 2 {R) must be real. Since the model of fi{R) 0 also 
satisfies this requirement during the BBN epoch, it looks 
like a possible cosmological model. 

The model fi{R) is, however, excluded since it can¬ 
not describe the late Universe, i.e., the ACDM model, 
whose energy density is dominated by the dark energy 
(A) and cold dark matter (CDM). Astronomical obser¬ 
vations indicate that the Ricci scalar becomes negative 
in the late Universe, which can be described with the 
CDM and dark-energy-dominated universe. For exam¬ 
ple, in the standard cosmological model, the Ricci scalar 
[Eq. (fT5]) [ is given by 


— R = (p — 3p) -I- 4A 

0 (radiation dominated epoch) 

= (matter dominated epoch) (23) 

4A (A dominated epoch) 


Since the present Universe is explained by a negative R 
value, the negative R should be consistently accommo¬ 
dated in the cosmological model. If the R value during 
BBN is positive as in the case of /i (i?), a transition from 
i? > 0 to i? < 0 must occur in the Universe between the 
BBN and the present epochs. While the f 2 {R) model 
can describe this late Universe, the fi{R) model cannot 


Qa — 


(6a)^“ 


— lOa^ -I- 8a — 1 
2|2a- 


(25) 


We show that for any fixed temperature T, the ex¬ 
pansion rate is larger for a larger value of a. For that 
purpose, we define a new parameter. 


h{a) 


/ H M4 
V Mp y Sttp 


4:06 

= or 9 


-1 

06 


/aN2“ 2|2a-l|^“^“ 
V6/ -10a2-h8a-l' 


(26) 


The derivative of the h{a) function is given for a y: 1/2 

by 


h'{a) 


/ON 2a 4|2a-l|^”^“ 

V6/ -10a2-F8a-l 


X 


In 


a 

12|a- 1/2| 


10 (a-2/5) 
-Wa^ + 8a-l 


■ (27) 


In the allowed parameter region of Eq. pl|l and 0.4347 < 
a < 1/2, the inequality h'{a) > 0 holds. The h{a) and 
the Hubble rate are, therefore, monotonically increasing 
functions of a in the vicinity of a = 1/2. 


III. BBN CALCULATION 

The public BBN calculation code ii is utilized and 
modified in this calculation. We updated rates of reac¬ 
tions related to nuclei with mass numbers <10 using the 
JINA REACLIB Database 0[ (the latest version taken 
in December 2014). The neutron lifetime is the cen¬ 
tral value of the Particle Data Group, 880.3 ± 1.1 s (6[. 
The baryon-to-photon ratio is (6.037± 0.077) x 10“^° [71, 
corresponding to the baryon density determined by the 
Planck observation of the cosmic microwave background, 
= 0.02205 ± 0.00028 [![. 
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In general, BBN codes take in only two independent 
equations from the equations of motion and the energy 
conservation equation. For example, the Kawano code 
Q uses equations of the Hubble expansion rate and the 
time derivative of temperature, i.e., dT / dt. In the present 
modified gravity model, the Hubble expansion rate is 
given by Eq. m- The energy conservation equation, 
i.e., Eq. (ini), is, on the other hand, not changed from that 
in the standard BBN (SBBN) model. We can, therefore, 
use the same equation for the time evolution of tempera¬ 
ture as that in the SBBN model [Eq. (D.26) in Ref. il- 
Then, only one modification of the Hubble rate to the 
code is required for BBN network calculations. 


IV. OBSERVED LIGHT ELEMENT 
ABUNDANCES 

Calculated BBN results are compared to the following 
observational constraints on light element abundances. 
The primordial ^He abundance is estimated with obser¬ 
vations of metal-poor extragalactic H H regions. We use 
the latest determination of Tp = 0.2551 ±0.0022 Q. The 
primordial D abundance is estimated with observations 
of metal-poor Lyman-a absorption systems in the fore¬ 
ground of quasistellar objects. We use the weighted mean 
value of D/H= (2.53 ± 0.04) x [l3|. ^He abun¬ 
dances are measured in Galactic H H regions through 
the 8.665 GHz hyperfine transition of ^He+. These are, 
however, not the primordial abundances but present val¬ 
ues which have also been also affected by Galactic chem¬ 
ical evolution, including production and destruction in 
stars. Nevertheless, it is very hard to reduce elemental 
abundances in the whole Galaxy, which have increased 
once by a significant factor since almost all of gas in the 
Galaxy needs to be incorporated into the stars multi¬ 
ple times and experience nuclear destruction reactions. 
We then adopt only the upper limit from the abundance 
^He/H=(1.9 ± 0.6) x 10“^ [ll| in Galactic H H regions. 
We note that this is neither a significant nor a direct limit 
on the primordial abundance and should be considered to 
be just a rough guide. The primordial ^Li abundance is 
estimated with observations of Galactic metal-poor stars. 
We use the abundance log(^Li/H)= —12 ± (2.199±0.086) 
derived in a 3D nonlocal thermal equilibrium model [l^ . 
®Li abundances in Galactic metal-poor stars have also 
been measured. We adopt the least-stringent 2 tr up¬ 
per limit of all limits for stars reported in [l^, i.e., 
®Li/H=(0.9±4.3) X 10-12 for the G64-12 (nonlocal ther¬ 
mal equilibrium model with five free parameters). 


V. RESULT 

Figure [T] shows the abundances of ^He (Yp; mass frac¬ 
tion), D, ^He, I’Li, and ®Li (number ratio relative to H) 
as a function of the power-law index of the scale factor 
a — 1/2 or the index n — 1 in the f{R) function. The left 


n-1 (xlO®) 




^ ^ ^ ^ T 
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- 


- 

- 


1 

1 
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a-1/2 (xlO®) 


FIG. 1: ^He mass fraction Yp and number abundance ratios of 
D, ®He, and ®Li relative to H as a function of the power- 
law index of the scale factor a —1/2 or the f{R) function n—1. 
Solid curves show calculated results for the f{R) oc i?" model. 
The solid smaller and larger boxes of "^He and D abundances 
correspond to the 2a and 4 (t limits, respectively, from adopted 
observational constraints. The dashed box corresponds to the 
2cr limits of ^Li abundance. The horizontal lines with down¬ 
ward arrows of ®He and ®Li abundances show observational 
upper limits. The result on the vertical line for a = 1/2 or 
n — I is for the SBBN model. Note that the f{R) function for 
a < 1/2 was corrected in this paper so that the cosmological 
model connects to the present ACDM model. 


half of the parameter region, i.e., a < 1/2, is consistent 
with the present ACDM model when the function f{R) 
[Eq. ((Ml) ] is adopted, while it is inconsistent when the 
function fi{R) [Eq. (fTOl) ] is adopted (see Sec. ITTl) . 

It is seen that the 2a observational limits of "‘He and D 
abundances are slightly inconsistent with the theoretical 
values for the SBBN, while their 4cr limits are consis¬ 
tent. In this figure, however, relatively small uncertain¬ 
ties coming from adopted nuclear reaction rates are not 
taken into account. Also, observational abundance deter¬ 
mination may include some unknown systematic error. In 
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addition, the disagreements of abundances are of the or¬ 
der of 10 % at most. The slight disagreements, therefore, 
do not seem so important at the moment, although they 
could be meaningful defects in the SBBN model in the 
future. In contrast to the disagreements, the ^Li abun¬ 
dance in the SBBN model is larger than observational 
limits by a factor of ^ 3 — 4. The SBBN values of ^He 
and ®Li abundances are, on the other hand, consistent 
with observational upper limits. 

As proven in Sec. m the cosmic expansion rate in the 
modified gravity model is increased when a increases. 
In particular, the expansion rate for a > 1/2 is always 
larger than in the SBBN model. All curves are under¬ 
stood as a result of a changed expansion rate as follows. 
First, when the expansion rate is larger, the freeze-out 
of weak reactions occurs earlier. The neutron abundance 
remaining after the freeze-out is then higher. Second, the 
time interval between the freeze-out and the ^He synthe¬ 
sis is shorter because of faster cosmic expansion. Neutron 
abundances are larger because of the above two reasons. 
Almost all neutrons are processed to form "^He nuclei at 
the ^He synthesis epoch. The ^He abundance is therefore 
larger for larger values of a. 

Because of the larger expansion rate, the reaction 
^H(n, 7 )^H also freezes out earlier. The relic neutron 
abundance right after the ^He synthesis then becomes 
higher. This higher neutron abundance significantly af¬ 
fects abundances of other light nuclei. D is predom¬ 
inantly produced via ^H(n, 7 )^H. The higher neutron 
abundance then leads to higher D abundance. is 
produced via and destroyed via n)^He. 

The enhanced D abundance leads to a higher abun¬ 
dance by a higher production rate. ^He is produced via 
n)^He and destroyed via ^He(n,p)^H. The some¬ 
what higher D abundance leads to a higher production 
rate, while the higher neutron abundance leads to a sig¬ 
nificantly higher destruction rate. Eventually, the ^He 
abundance is slightly higher. The primordial abun¬ 
dance is the sum of and ^He produced during the 
BBN. Long after the BBN, ^He nuclei /3-decay into 
nuclei. The final abundance of ^He is larger than that of 
by about 2 orders of magnitude in SBBN. Therefore, 
the primordial abundance predominantly reflects the 
larger ^He abundance during BBN. 

^Li is prodnced via ^He(t, 7 )^Li and destroyed via 
^Li(p, a)^He. The ^Li abundance is then higher be¬ 
cause of the higher T abundance. ^Be is produced 
via "^He(^He, 7 )^Be and destroyed via ^Be(n,p)^Li. A 
slightly higher abundance of ^He and a considerably 
higher abundance of the neutron leads to a smaller ^Be 
abundance. The primordial ^Li abundance is the sum 
of ^Li and ^Be produced during the BBN. Long after 
the BBN, ^Be nuclei recombine with electrons and are 
transformed to ^Li nnclei via the electron capture pro¬ 
cess. The abundance of ^Be is larger than that of ^Li in 
SBBN. Therefore, the larger expansion rate results in a 
smaller primordial ^Li abnndance. ®Li is produced via 
"^He(d, 7 )®Li and destroyed via ®Li(p, Q;)^He. The higher 


D abundance leads to a higher ®Li abundance. 

We note that trends of theoretical curves in the mod¬ 
ified gravity model are somewhat similar to those in the 
model including additional components of energy density, 
e.g., sterile neutrinos and a primordial magnetic field, 
within the framework of Einstein’s general relativity. A 
more illnstrative and detailed explanation can be found 
for effects of the magnetic field on all light element abun¬ 
dances in Ref. M- 

When the 4 (t limit of ^He is used, the constraints are 
derived, 

- I X lO"*^ < (a - 1/2) < 5 X lO"'^ 

-2xI0"'^< (n-I) < 10“^ (28) 

When the 4cr limit of D is used, however, more stringent 
constraints are derived, 

- 3 X IQ-'^ < (a - 1/2) < 2 X I0“® 

-exIQ-'^^ (n-I) <4xI0“®. (29) 

We find that the cases of a > 1/2 and a < 1/2 are con¬ 
strained most stringently from the abundance constraints 
on D and ^He, respectively, for the following reasons. 
The 4(7 range of D abundance is barely consistent with 
the SBBN result, and the upper limit is very close to 
the SBBN value. Since the D abundance increases with 
increasing a ot n, the theoretical curve easily deviates 
from the observational limit for a > 1/2 (Fig. [T|). On 
the other hand, the lower limit of ''’He abundance is very 
close to the SBBN value. Since the ''He abundance also 
increases with increasing a or n, the theoretical curve 
easily deviates from the observational limit for a < 1/2. 

Additionally, we find that the abundances of ^He and 
6Li are not sensitive to the change of a, and that they 
are within the adopted limits in the whole parameter 
region shown in Fig. [TJ Althongh the ^Li abundance 
approaches the observed abundance level with increasing 
a, it is always ontside of the adopted limits in Fig. [TJ 

Fignre 121 shows the likelihood function for the power- 
law index a—1/2 or n—1 (solid curve). In this estimation, 
the likelihood function is defined as 


L(a) = LD(a)L4jjg(a), (30) 


where the likelihood functions for respective nuclear 
abundances (i =D and ^He) are given by 


Li(a) 




obs 



[^,th(o^) ^,obs] 


2(7 


2 

2 ,obs 


(31) 


where Li,th(cr) is the theoretically calculated abundance 
of i, and Li,obs and cr^ obs are the central value and 1 a 
error, respectively, of the adopted observational abun¬ 
dance of i. The solid curve has been normalized as 
/ L{a)da = 1. 

From the combination of the constraints on D and ''He 
abundances, we derive a 95 % C.L. on the parameters as 

{a-1/2) = (-0.43 ± 0.59) X 10"'' 

(n-1) = (-0.86 ± 1.19) X 10"''. (32) 
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n-l (xlO^) 

-3-2-10 1 2 



FIG. 2: Likelihood function as a function of the power-law 
index a — 1/2 or n — 1 (solid curve). This curve is estimated 
using the observational constraints on the D and ^He abun¬ 
dances. The parameter region bounded by vertical dashed 
lines is the 95 % C.L. The vertical solid line on a = 1/2 or 
n = 1 corresponds to the SBBN model. 


This parameter region is enclosed by vertical dashed 
lines. The solid vertical line corresponds to the SBBN 
model described by the general relativity {a = 1/2 or 
n = 1). It is found that the deviation of the parame¬ 
ter a from 1/2 is constrained to be less than 0(10“^). 
This constraint also indicates that even in the generalized 
model of the f{R) gravity, the case of u = 1 is very likely 
based upon the comparison of theoretical and observa¬ 
tional nuclear abundances. If only the region of a > 1/2 


is allowed, as in the case of the f 2 {R) function [Eq. (fTTl) ]. 
the amplitude of a—Xj^ is constrained relatively strongly. 


VI. SUMMARY 

In this study we revisited effects of a modified gravity 
on BBN. The model is based on the /(i?) oc i?" term in 
the action and the assumption of the scaling for the time 
evolution of scale factor a{t) oc in a homogeneous and 
isotropic universe. The f{R) functions were constructed 
so that they describe the accelerated expansion of the 
present Universe successfully. We utilized a nuclear reac¬ 
tion network code and calculated all light element abun¬ 
dances in the modified gravity model. We compared the 
calculations with astronomical observations of primor¬ 
dial elemental abundances and found that the param¬ 
eters are constrained to be (a —1/2) < 0(10“''’) and 
(u — 1) < 0(10“^) mainly from the limits on primordial 
D and ''He abundances. 


Acknowledgments 

This work was supported in part by the Na¬ 
tional Research Foundation of Korea (NRF) (Grants 
No. NRF-2012R1A1A2041974 and No. NRF- 
2014R1A2A2A05003548). S.K. was supported by the 
Basic Science Research Program through the NRF 
funded by the Ministry of Education (No. NRF- 
2014R1A1A2059080). 


[1] G. Lambiase and G. Scarpetta, Phys. Rev. D 74, 087504 
(2006). 

[2] J. U. Kang and G. Panotopoulos, Phys. Lett. B 677, 6 
(2009). 

[3] L. Kawano, NASA STI/Recon Technical Report N 92, 
25163 (1992). 

[4] M. S. Smith, L. H. Kawano, and R. A. Malaney, Astro- 
phys. J. Suppl. 85, 219 (1993). 

[5] R. H. Cyburt et al, Astrophys. J. Suppl. Ser. 189, 240 

( 2010 ). 

[6] K. A. Olive et al. (Particle Data Group Gollaboration), 
Chin. Phys. C 38, 090001 (2014). 

[71 H. Ishida, M. Kusakabe, and H. Okada, Phys. Rev. D 90, 
083519 (2014). 

[8] P. A. R. Ade et al. (Planck Collaboration), Astron. As¬ 
trophys. 571, A16 (2014). 

[9] Y. 1. Izotov, T. X. Thuan, and N. G. Guseva, Mon. Not. 
Roy. Astron. Soc. 445, 778 (2014). 

[10] R. Cooke, M. Pettini, R. A. Jorgenson, M. T. Murphy, 
and C. C. Steidel, Astrophys. J. 781, 31 (2014). 

[11] T. M. Bania, R. T. Rood, and D. S. Balser, Nature 415, 


54 (2002). 

[12] L. Sbordone et al., Astron. Astrophys. 522, A26 (2010). 

[13] K. Lind, J. Melendez, M. Asplund, R. Collet, and 
Z. Magic, Astron. Astrophys. 544, A96 (2013). 

[14] M. Kawasaki and M. Kusakabe, Phys. Rev. D 86, 063003 

( 2012 ). 

[15] When another choice of the signs in the metric, goo = 
— 1 and gij = a{t)^Sij, is adopted, the formulation is 
somewhat different from the present case. The equations 
of motion are different from Eqs. © and ®, and the 
Ricci scalar, denoted as R*, is opposite in sign to R [Eq. 
(nsii. We find that, in general, the different choices of 
the sign in the metric result in different equations for the 
evolution of p and p [cf. Eqs. (1131) and (1141) 1 if the same 
function, i.e., f{R) and f{R*), respectively, is used. The 
functions /i(R*) = -{-R*/AY and fiiR*) = (R’/A)" 
then correspond to fi{R) and / 2 (R), respectively, in the 
present metric case. In the case of n = 1, however, it 
follows that fi{R) = f 2 {R) = R/A, and we obtain the 
same equations for p and p from Eq. ([2} independently 
of the choice of metric. 






